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ABSTRACT: In the context of the decoherent histories approach to quantum theory, it 
is shown that a class of macroscopic configurations consisting of histories of local densities 
(number, momentum, energy) exhibit negligible interference. This follows from the close 
connection of the local densities with the corresponding exactly conserved (and so exactly 
decoherent) quantities, and also from the observation that the eigenstates of local densities 
(averaged over a sufficiently large volume) remain approximate eigenstates under time 
evolution. The result is relevant to the derivation of hydrodynamic equations using the 
decoherent histories approach. 



The key step in many studies of emergent classicality from quantum mechanics is 
the demonstration that certain types of quantum states of the system exhibit neghgible 
interference. Initial superpositions of such states may therefore be effectively replaced 
by statistical mixtures. This, loosely speaking, is decoherence, and has principally been 
demonstrated for the situation in which there is a distinguished system, such as a particle, 
coupled to its surrounding environment [1]. 

Most generally, decoherence typically comes about when the variables describing the 
entire system of interest naturally separate into "slow" and "fast", whether or not this 
separation corresponds to, respectively, system and environment^. If the system consists 
of a large collection of interacting identical particles, as in a fluid for example, the natural 
set of slow variables are the local densities: energy, momentum, number, charge etc. These 
variables, in fact, are also the variables which provide the most complete description of 
the classical state of a fluid at a macroscopic level. The most general demonstration of 
emergent classicality therefore consists of showing that, for a large collection of interacting 
particles described microscopially by quantum theory, the local densities become effectively 
classical. Although one might argue that the system-environment mechanism might play a 
role, since the collection of particles are coupled to each other, decoherence comes about in 
these situations for a different reason: it is because the local densities are almost conserved 
if averaged over a sufficiently large volume [3] . Hence, the approximate non-interference of 
local densities is due to the fact that they are close to a set of exactly conserved quantities, 
and exactly conserved quantities obey superselection rules. 

Intuitively appealing though this argument is, it is clearly a quantitative issue. The 
object of this letter is to show that, under certain reasonable conditions, local densities 

t See Ref.[2] for a discussion of the conditions under which the total Hilbert space may be 
written as a tensor product of system and environment Hilbert spaces 
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averaged over a sufficiently large volume are indeed approximately decoherent as a result 
of their close connection to exact conservation. 

We will approach the question using the decoherent histories approach to quantum 
theory [3-5], which has proved particularly useful for discussing emergent classicality in a 
variety of contexts*. The central object of interest is the decoherence functional, 



The histories are characterized by the initial state and by the time-ordered strings of 
projection operators Cq = Pa^l^n) • • • Paiih) (where a denotes the string of alternatives 
CKl, Q!2 • • -CKn)- Intuitively, the decoherence functional is a measure of the interference be- 
tween pairs of histories a, a'. When it is zero for a a', we say that the histories are 
decoherent and probabilities p{a) = D{a,a) obeying the usual probability sum rules may 
be assigned to them. Although not addressed here, one can then ask whether these prob- 
abilities are strongly peaked about trajectories obeying classical equations of motion. For 
the local densities, these equations will be hydrodynamic equations, and these and closely 
related aspects of emergent classicality have been pursued at greater length elsewhere 



We consider the class of systems which are described at the microscopic level by a 
Hamiltonian of the form 



For definiteness, we will concentrate on the case of a dilute gas with short-range inter- 
actions, but it will be clear that the physical ideas are reasonably general. The local 

* The extent to which the approach fully explains emergent classicality has been criticized 
[6]. This paper concerns the mathematical properties of the approach, as it currently 
stands, and adds nothing to that debate. See Ref.[7] and references therein for further 
discussion. 




(1) 



Refs.[2,7-9]. 
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densities of interest are the number density n(x), the momentum density g(x) and the 
energy density /i(x), defined by, 

n(x) = ^ 5(x-q,-) (3) 
3 

g(x) = ^ Pj 5(x-q_^) (4) 
3 

i \ i>3 I 
(suitably ordered, in the quantum case). We are interested in local densities smeared over 

a volume V . The eff'ect of this is to replace the delta functions with a window function, 

denoted which is zero outside V and 1 inside. It is also useful to work with the Fourier 

transforms of the local densities, denoted n(k), 5r(k), /?,(k). So, for example, the local 

number density at wavelength k is 

n(k) = J] (P^-"^ (6) 

3 

Exact conservation is obtained in the limit k = |k| — » 0, or V oo in (3)-(5). 

We would like to compute the decoherence functional for histories consisting of projec- 
tions onto the operators (3)-(5). (The construction of the projectors is described in more 
detail in Ref.[7]). In the case of exact conservation. A; = 0, we have exact decoherence 
simply because the projectors in Eq.(l) all commute with H and with each other [10]. 
Our main task is therefore to show that as k increases from zero there is still a non-trivial 
regime in which decoherence is approximately maintained. A significant result of this 
type has been established already by Calzetta and Hu for the case of local temperature 
T{x) obeying the diffusion equation [9]. They took their initial state to be close to the 
equilibrium state, whereas here, by contrast, initial macroscopic superposition states are 
considered. 
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We begin by rewriting the exact conservation case in a simple way that makes its 
generalization to locally conserved quantities more apparent. Suppose the histories are 
projections onto some conserved quantity, Q. Let the initial state be a superposition of 
eigenstates of Q, 

m = ^{\a) + \b)) (7) 

where {a\b) = and 

Q\a)^a\a), Q\b) = b\b) (8) 

Since the P^s are projections onto Q, Pa either annihilates or preserves |a) and \b). Take 
the case of a history with just two moments of time (the generalization to more times is 
trivial) . The only non-zero off-diagonal terms of the decoherence functional are of the form 

D{a,a!) = ^Tr (^Pa,e-i^*\a){b\ei^*^ 

= ^TT{Pa,\at)m (9) 

But Q is conserved, hence [Pa2 , H] = and 

PaM = Pa2e~i^*\a) 

= e-i^^Pa,\a)^\at) (10) 

(or equals zero if a2 does not correspond to a). It follows that 

D{a, a!) = ^Tr {Pa,\at){bt\) 

= {bt\at) = {b\a) = (11) 

and therefore we have decoherence. 

Now suppose that the operator Q is one of the local densities (3)-(5), so is no longer 
exactly conserved. The steps up to Eq.(9) still hold. But to go further, we need to 
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know how the eigenstates of the local densities behave under time evolution. A reasonable 
supposition, which will be justified, is the following. Let us suppose that under time 
evolution, the eigenstates of Q remain approximate eigenstates. That is, we initially have 



i.e., the state remains strongly peaked in the variable Q under time evolution. The states 
are then approximate eigenstates of the projectors, so that in place of Eq.(lO), we have the 
approximate result, P^gl^t) ~ Wt) (or equals zero) as long as the width of the projection 
is much greater than the uncertainty (AQ)^. Hence Eq.(ll) follows approximately, and 
we get approximate decoherence to the extent that the approximation (12) holds. 

The key point is therefore the following: approximate decoherence is assured for histo- 
ries of operators Q whose eigenstates have the property that they remain strongly peaked 
in Q under time evolution, as characterized by (12). To demonstrate decoherence of the 
local densities, therefore, we need only find their eigenstates, and show that they satisfy 
the localization property (12) under time evolution. (Note, incidently, that the above argu- 
ment actually assures decoherence of any variables Q satisfying the localization property. 
The particular significance of the local densities is that they are continuous functions of the 
coarse graining scale k, so are guaranteed to satisfy the requisite property if k is sufficiently 
close to zero.) 

Since the three operators (3)-(5) do not commute, exact simultaneous eigenstates 
cannot be found. However, a useful class of approximate eigenstates of all three operators 
are the states consisting of products of N identical terms. 



(7), but under evolution to time t, Q\at) ~ {Q)\at) or, more precisely. 



(AQ)' 



« 1 



(12) 




(13) 
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These may be shown to be eigenstates (of the local number density, for example) by 
observing that the object (An(x))^/(?2(x))^ goes like 1/N for large (see Ref.[7], for 
example). It is essentially the central limit theorem (see also Ref.[ll]). For the number 
and momentum density it relies on the fact that they are sums of identical one-particle 
operators. For the local energy density, it additionally requires the smearing volume to 
be sufficiently large, compared to some lengthscale indicated by the interactions. Some 
tuning of the state jV') can be carried out to ensure that (13) is an optimal approximate 
eigenstate of all the local densities but this will not be done here. (Also, the passage to 
exact eigenstates of n(k), ^(k), /j,(k) as A; ^ can be seen explicitly if the one-particle 
states are taken to be one-particle momentum eigenstates). 

The question is now what happens to the eigenstates (13) of the local densities under 
time evolution by the Hamiltonian (2). Consider first the trivial but enlightening case 
in which there no interactions. In this case, the time evolved eigenstates \at) remain of 
the product form (13), so they are still approximate eigenstates of the local densities (but 
with a time-evolved eigenvalue) for the same reasons as above. Hence there is approximate 
decoherence. 

Decoherence in the non-interacting case comes about for two reasons. First, it is due 
to the fact that a state of the form (13) will remain strongly peaked about the average 
values of the local densities, n(x), 5r(x), /i(x) under time evolution, and thus the state is 
essentially undisturbed by the projectors (as long as their widths are sufficiently large). 
Secondly, it is due to the almost trivial fact that the orthogonality of the two elements of 
the initial state is preserved by unitary evolution. This second fact is important because 
the first one is not always sufficient to guarantee decoherence. Although the state remains 
strongly peaked about the average values of the local densities, these average values do not 
necessarily obey deterministic equations. In the case of histories characterized by number 
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density only, for example, (?^(x)) at time t is not uniquely determined by (n(x)) at the 
initial time (in the state (13)). That is, in Eq.(9), \at) and \bt) may in fact be peaked 
about the same value of number density, even though the initial values are different. The 
decoherence is therefore not in fact due to an approximate determinism (such as that 
used in the phase space histories of Omnes [5]). It is necessary only that the evolved 
states are essentially undisturbed by the projectors, and therefore that the two orthogonal 
components of the initial state are eventually overlapped at the final time, as in Eq.(ll), 
to give zero. 

The next and most important task is to show that the above story is in fact still 
true, with qualifications, in the presence of interactions. The complete description of 
N interacting particles is complicated but we can make some progress by making two 
assumptions which are standard in kinetic theory [12]. It is notationally convenient in what 
follows to work with a Wigner function, rather than quantum state. Hence associated with 
the full A'"-particle wave function is an AT-particle Wigner function W7v(pi, Ql, • • • PiV, QAr)- 
Our first assumption is that the three-particle correlations are negligible. It means that all 
the physics is contained in the one and two-particle reduced Wigner functions, Wi(pi, qi) 
and W2(pi,qi,P2,q2)- 

We again take as our initial state the approximate eigenstate (13), and let it evolve, so 
correlations will develop. The degree to which the particles become correlated is contained 
in the two-particle distribution W2 of the evolved eigenstate. On general grounds, we 
expect that the inter-particle correlations will only be important on some length scale L, 
and beyond that length scale, they will be uncorrelated. That is, we will assume that 

w^2(pi, qi, P2, q2) ~ w^i(pi, qi)Wi(P2, q2) (14) 

for |q2 — qil > and otherwise W2 will have a form indicating non-trivial correlations. 
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This is our second assumption. It is physically reasonable for uncorrelated initial states 
of the form (13) with a short range interaction (and it is in fact a key assumption in the 
derivation of the Boltzmann equation [12]). It would not of course be an appropriate as- 
sumption for correlated initial states, but the point is that we are interested in approximate 
eigenstates of the local densities, and a useful class of such states have the uncorrelated 
from (13). 

Given the above assumptions, it is now reasonably straighforward to argue that the 
state is still strongly peaked about the average values of the local densities, as long as 
V » L^. For example, for the number density, we have 

(n(x)) = (q,- -^)) = N I d\ p(q) (15) 

. JV 

where p(q) is the one-particle probability distribution of q (obtained by integrating the 
one-particle Wigner function over p). Similarly, 

(n2(x)) = (^^(^i - ^)^v{^ - x)) 

= N{5v) + (iV^ - iV)(5y(qi - a;)5y (q2 - x)) (16) 

where we have used 5y = and also an assumption of identical particles to reduce the 
sum over j, I to particles labeled 1 and 2. We now have 

(An(x))2 =(n2(x)) - (n(x))2 

=iV2 ((5^(qi - x)5y(q2 - x)) - {Sy?) 

+ N {{5y) - (V(qi - x)5v.(q2 - x))) (17) 

If there is no correlation at all between the particles, the coefficient of N'^ would vanish, 
so (An(x))^/(n(x))^ would go like which goes to zero as A?" — > oo. This is the 

standard central limit theorem result indicated earlier for the non-interacting case. With 
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interactions, the coefficient of N'^ is no longer zero. We now need to show, therefore, 
that this term is still sufficiently small for (Ari(x))^/(n(x))^ to remain small as AT — > oo. 
Introducing the two-particle distribution p(qi,q2) (obtained by integrating pi,P2 out of 
W2), it is readily shown that the leftover terms as AT —> 00 are 

(An(x))2 ^ Jy d^qi Jy d^q2 (p(qi, q2) -p(qi)p(q2)) .^gs 

This is clearly zero if there are no correlations. In the interacting case we use the assump- 
tion (14), which implies that 

p(qi, q2) ~ p(qi)p(q2) (i9) 

for |qi — q2| > -^5 and otherwise non-trivial correlations exist. Hence the integral in the 
numerator takes contributions only from the region |qi — q2| < L. 

To see that (18) is small, note that in the numerator, the integral is over a volume 
in the six-dimensional two particle configuration space. If F << L^, the factorization of 
p(qi) q2) for |qi — q2| > L makes no difierence, since qi and q2 can never be far enough 
apart in the integrand (assuming V is regular in shape). However, if V » L^, the F^- 
sized integration region is substantially reduced in size to F x L^. On dimensional grounds 
the numerator is therefore proportional to a number of order VL^, and the denominator 
to (perhaps with other factors common to both). This means that 

, / ~ -77 (20) 

This order of magnitude estimate becomes exact if we assume that the probabilities are 
constant in the region of non-trivial correlation (another common assumption of kinetic 
theory [12]). Hence the state will be strongly peaked about the average of n(x) if F >> L^. 

It is possible to see on physical grounds why one expects a result of the form (20) to hold 
quite generally. In the non-interacting case we used the central limit theorem result that 
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(An) / (n) goes like 1 /N . In the interacting case, the state is no longer of the product form 
(13), but an analagous result still holds. The point is that the correlations that develop 
extend only over a (typically small) volume of size L^, so the system breaks up into a large 
number of essentially identical uncorrelated regions of this size. Therefore each smearing 
volume V , if much greater than L^, contains of order V/L^ identical uncorrelated regions 
each of which contribute equally to the local density averaged over V . Loosely speaking, 
a central limit theorem-type result again applies, not to the uncorrelated particles in 
the same state, but to the V/L^ uncorrelated regions. So 1/A^ is replaced by /V in the 
central limit theorem, and hence the above result. 

Similar results hold for the local momentum and energy density. We have therefore 
demonstrated the desired result: a class of eigenstates of the coarse-grained local densities 
remain approximate eigenstates under time evolution as long as the smearing volume is 
much greater than the correlation volume of these states. Decoherence of these variables 
then follows. More details of this work, including a discussion of the approach to local 
equilibrium and the emergence of hydro dynamic evolution equations, will be published 
elsewhere [13]. 
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